raw=OMFIT['ASIPP_IM']['EASTdatas']['Raw_datas']
shot_list=range(142616,142647,1)
time_list=[3.5,5,6,8]
for shot in shot_list:
    if str(shot) not in raw:
        raw[str(shot)]=OMFITtree("")
    for time in time_list:
        if str(time) not in raw[str(shot)]:
            raw[str(shot)][str(time)]=OMFITtree("")
        # load gfile
        gfile_path='/gpfs/scratch/hyueheng/mgfile/'+str(shot)+'/g'+str(shot)+'.00'+str(int(time*1000))
        if os.path.isfile(gfile_path):
            raw[str(shot)][str(time)]['mgfile']=OMFITgeqdsk(gfile_path)
        # load ti
        Tifile_path='/gpfs/scratch/hyueheng/Ti/'+'/Ti_'+str(shot)+'_'+str(time)+'.mat'
        #print gfile_path
        if os.path.isfile(Tifile_path):
            raw[str(shot)][str(time)]['Tifile']=OMFITmatlab(Tifile_path)
        # load te
        Tefile_path='/gpfs/scratch/hyueheng/ECE/'+'/Te_ECE_'+str(shot)+'_'+str(time)+'.mat'
        #print gfile_path
        if os.path.isfile(Tefile_path):
            raw[str(shot)][str(time)]['Tefile']=OMFITmatlab(Tefile_path)
            

